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Abstract 

We explore the cosmological solutions of a recently proposed extension of 
General Relativity with a Lorentz-invariant mass term. We show that the 
same constraint that removes the Boulware-Deser ghost in this theory also 
prohibits the existence of homogeneous and isotropic cosmological solutions. 
Nevertheless, within domains of the size of inverse graviton mass we find 
approximately homogeneous and isotropic solutions that can well describe the 
past and present of the Universe. At energy densities above a certain crossover 
value, these solutions approximate the standard FRW evolution with great 
accuracy. As the Universe evolves and density drops below the crossover value 
the inhomogeneities become more and more pronounced. In the low density 
regime each domain of the size of the inverse graviton mass has essentially non- 
FRW cosmology. This scenario imposes an upper bound on the graviton mass, 
which we roughly estimate to be an order of magnitude below the present-day 
value of the Hubble parameter. The bound becomes especially restrictive if 
one utilizes an exact self- accelerated solution that this theory offers. 

Although the above are robust predictions of massive gravity with an ex- 
plicit mass term, we point out that if the mass parameter emerges from some 
additional scalar field condensation, the constraint no longer forbids the homo- 
geneous and isotropic cosmologies. In the latter case, there will exist an extra 
light scalar field at cosmological scales, which is screened by the Vainshtein 
mechanism at shorter distances. 



1 Introduction and Summary 



The purpose of this work is to study the cosmology of General Relativity (GR) with 
an explicit Lorentz invariant mass term (massive gravity or massive GR). Histori- 
cally, it has been difficult to construct a nonlinear theory of massive gravity that 
would describe no more than the five degrees of freedom, required for the massive 
spin-2 state by the representations of the Poincare group. This situation was re- 
cently transformed by the proposal in Ref. [1] of a theory of massive gravity with 5 
degrees of freedom. 

This theory was shown to be free of the sixth degree of freedom (the so-called 
Boulware-Deser (BD) ghost [2]) to all orders in the decoupling limit (DL) in Refs. 
[1, 3], where it was also shown to be ghost-free away from the DL up to and including 
quartic order in nonlinearities [1]. These arguments were recently generalized in [4] 
to a complete nonlinear proof of the absence of ghosts, away from the decoupling 
limit, in the ADM/Hamiltonian formalism. The same result can be reached in both 
the Stiickelberg (see Ref. [5]) and helicity formalisms (see Ref. [6]). In each language, 
there exists a constraint that eliminates one degree of freedom which otherwise would 
have been the BD ghost. 

In this work we will show that the very same constraint that eliminates the BD 
ghost in massive gravity [1], also forbids homogeneous and isotropic cosmological 
solutions (FRW cosmologies). For solutions with FRW symmetries, the all-orders 
constraint can straightforwardly be seen to prohibit any time evolution, leaving 
Minkowski space as the only vacuum solution which is consistent with homogeneity 
and isotropy. 

This result raises the question: How could the non-FRW cosmologies of massive 
gravity recover the FRW solutions of GR in the massless limit? For this one should 
rely upon the Vainshtein mechanism by which massive GR is expected to recover 
GR in the m — > limit [7]. Although the original theory in which the Vainshtein 
mechanism was proposed contains a BD ghost, the mechanism itself seems to be 
universal, and has been established in other models where the BD ghost is not 
present [8, 9, 10, 11]. Moreover, it was shown that the mechanism is operative for 
spherically symmetric solutions in the massive GR theories discussed here, at least 
for a certain choice of the three a priori free parameters of the theory (i.e., the 
graviton mass and two arbitrary constants), [12, 13]. 

Assuming that the Vainshtein mechanism is at work, one would expect to find 
in massive GR cosmological solutions that are more and more homogeneous and 
isotropic as the value of the graviton mass is taken to zero. If this is the case, then 
the fact that massive gravity leads to non-FRW solutions will not immediately rule 
it out via observations, but rather just place a constraint on the magnitude of the 
mass of the graviton, to be consistent with known constraints on homogeneity and 
isotropy. 

To see a close connection between the Vainshtein mechanism and cosmology, 
consider matter of constant density p stored in a sphere of radius R. The Vainshtein 



1 



radius of such a source is 




where r g = 2MGn is the gravitational radius of the source of mass M = ^rrR 3 p, 
and G N = (SirM^)' 1 . Furthermore, it is useful to introduce a notion of a crossover 
energy density, 

Pco = 3M^m 2 . (2) 

From (1) we conclude that any source with density above the crossover value (2), 
is characterized by the Vainshtein radius that is greater than the size of the source 
itself. For such sources gravity is close to that described by GR at distance scales 
<C r* (the Vainshtein regime), however, deviates significantly from GR at distance 
scales > r* (the vDVZ regime) [14] . Furthermore, sources with density below the 
crossover p co , are always in the vDVZ regime and their gravity differs significantly 
from GR. 

Let us now apply these observations to cosmology. Suppose we took a snap- 
shot of a universe at a certain stage of its evolution when the matter in it had an 
average energy density p (averaged, say, at scales greater than the Hubble scale 
H -1 at that epoch; here we suppose that inflation, or an alternative early universe 
framework, prepared such a state). Let us look into a 1/m-size domain. An ar- 
bitrary Hubble patch in this domain (i.e., a patch enclosed by a sphere of radius 
H^ 1 = (p/SMp^^ 1 / 2 ), that is far enough from the edges of the domain, is well- 
within the Vainshtein regime as long as p > p co . Then, cosmology within such 
Hubble patches can be approximated by the standard FRW metric of GR with 
small corrections. Hence, for p ^> p co , and well-within each 1/m-size domain, the 
early universe in massive GR would evolve as it does in GR, with some small cor- 
rections that vanish in the m — > limit. Such an expansion in each of these Hubble 
patches will last until the size of the patch, H' 1 , approaches the scale ~ 1/m, or 
equivalently, until p dilutes down to density of the order of p co . At scales larger than 
1/m gravitational interactions are expected to be screened. 

The same arguments should hold for the radiation dominated epoch, in which 
case p ~ T 4 (T being temperature) should be compared with p co . 

Requiring that the graviton mass be less then the Hubble parameter today, 
m < H , we find that p co < p c , where p c , is the present-day value of the critical 
density in the Universe. If so, then according to the above described scenario, the 
cosmological evolution of the early Universe (p ^> p co ) within each 1/m-size domain 
will mimic the FRW expansion with some accuracy. However, this will change 
significantly at densities p ~ p co . As long as the graviton mass m is sufficiently 
small, the observational tests of such cosmologies of massive GR would impose 
an upper bound on m. We estimate this bound to be approximately an order of 
magnitude smaller than Ho ~ 1CT 33 eV. 
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In general, the mass term introduces an effective stress-tensor in the Einstein 
equation. The backreaction of this term should be negligible in the Vainshtein 
regime, becoming dominant in the vDVZ regime. While generically we expect the 
above to be the case, interestingly enough, we find one particular exact solution for 
which the backreaction is described by a perfect fluid with the equation of state of 
dark energy, and the magnitude of the energy density/pressure set by m 2 M^, v This 
behavior is similar to the self-accelerated solutions of massive GR first found in Ref. 
[15] in the DL or its extension in Ref. [16], and to exact self-accelerated solutions 
obtained in Refs. [12, 17, 18]. 

If one utilizes this particular solution, then by the end of the Vainshtein regime 
the universe may become dominated by the self-accelerated solution. However, the 
latter is not regular at spatial infinity, and could only exist as a transient solution 
in space and time, matched upon the low-density inhomogeneous solution at larger 
scales 1 . Whether such a matching is possible, is not shown here; in principle the 
evolution could just bypass this solution and transition directly to a low density 
regime. Putting the question about the matching aside, however, the bound on 
the graviton mass becomes especially restrictive if the expansion of the universe is 
described by the self- accelerated solution (see discussions in section 3). 

In this paper we consider the cosmological evolution both in the Vainshtein and 
vDVZ regimes. In the Vainshtein regime, the metric to which the matter couples 
is homogeneous and isotropic with some small corrections, but it's the Stiickelberg 
sector that carries all the inhomogeneities. Moreover, in this regime, and for the 
self-accelerated solution, we will show the existence of the backreaction of the mass 
term that is small and mimics dark energy. What is not shown is that there is 
a matching between the Vainshtein (with or without self-acceleration) and vDVZ 
regimes 2 . 

The theory of Ref. [1] is the only potentially viable classical theory of Lorentz 
invariant massive GR with 5 helicity states. Its cosmology is unusual, and this paper 
is a first attempt at unfolding peculiarities of such a theory in a cosmological setup. 
Therefore, the majority of this paper is qualitative in character, where we emphasize 
certain universal aspects and set up a general framework in which such cosmologies 
can further be studied in details. A number of particular exact cosmological solutions 
are discussed in the appendices. 

1 Note that the solutions of Refs. [12, 17] exhibit singularities at finite values of the coordinates, 
and hence, should be matched to other solutions before reaching those points. 

2 To address this issue, one could consider a possibility that the matter/radiation that is being 
expelled from the bulk of the 1/ra-size domains, which are densely packed and adjacent to each 
other, gets accumulated near the boundaries of the domains. If the domains are well separated, 
or there is only one domain, density near the edge will be suppressed due to screening of gravity 
and free streaming. Different sccnaria arc determined by different initial conditions and need more 
detailed studies. 
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2 Massive GR and Cosmology 



For massive GR the action is a functional of the metric g^x), and four spurious 
scalar fields <f) a (x), a = 0,1,2,3; the latter are introduced to give a manifestly 
diffeomorphism invariant description [19, 20]. One defines a covariant tensor H^ v as 
follows: 

SV = d^d^rjab + , (3) 

where r] ab = diag(— 1, 1, 1, 1). The first term on the r.h.s. is nothing but the 
Minkowski metric in the coordinate system defined by a 's. Hence, gravity in this 
formulation is described by the tensor H^ u propagating on Minkowski space. In 
the unitary gauge all the four scalars <f) a (x) are frozen and equal to the correspond- 
ing space-time coordinates, (f> a (x) = x^5^. However, often it is helpful to use a 
non-unitary gauge in which a (x)'s are allowed to fluctuate. 

A covariant Lagrangian density for massive GR can be written as follows, 

C=^y^(R-^U(g,H)) , (4) 

where U includes the mass, and non- derivative interaction terms for H^ u and g^ u . 

A necessary condition for the theory to be ghost free in the DL is that the 
potential ^/^g ^(fiS H) be a total derivative upon the field substitution = g^ u — 
f]ixv = 0, <fi a = S^x^ — rj^d^Tr [3]. With this substitution, the potential yf^g U(g, H) 
becomes a function of = d^d^n and its various contractions. 

For instance, the following expression composed of is a total derivative 

4 2 2(H) ee [H] 2 -[H 2 ], (5) 

where we use the notations [n] = trll£, [n] 2 = (trn^) 2 , while [n 2 ] = trn^n^. 

Then, as argued in [1], the Lagrangian for massive GR that is automatically 
ghost free to all orders in the DL is obtained by replacing the matrix elements 
in the total derivative term (5) by the matrix elements of a tensor )C£, defined as 
follows: 

ty(g,H) = ^-VWV. (6) 

Here, the indices on K, should be lowered and raised by g^ v and its inverse respec- 
tively. This procedure defines the mass term (along with the interaction potential 
in the Lagrangian density) in massive GR 

£ = ^V^g [R - m> (wic; - ra 2 )] . (7) 

The matter and other fields are coupled to g^ as in GR. The above expression 
has no free parameters once the graviton mass is fixed. In general, however, there 
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exist other polynomial terms in /C with similar properties. These terms can be 
constructed straightforwardly by using the procedure outlined in Ref. [1]. In any 
dimensions there are only a finite number of total derivative combinations, made of 
II [11]. They are all captured by the recurrence relation [3]: 

An) _ V-V ( n ~ jO; [ n ml r (n-m) /gx 

^ cr --2^(-i) P 14* , (8) 

m=l v ' 

with = 1 and C^ T = [II]. This also guarantees that the sequence terminates, 
i.e., £^ = 0, for any n > 5 in four dimensions. The list of all nonzero total 
derivative terms starting with the quadratic one reads as, 

4 2 2(n) = [n] 2 -[n 2 ], (9) 
££(n) = [n] 3 -3[n][n 2 ] + 2[n 3 ], (10) 
4 4 2(n) = [n] 4 -6[n 2 ][n] 2 + 8[n 3 ][n] + 3[n 2 ] 2 -6[n 4 ]. (11) 

One can use the method of Ref. [1] to obtain the two other polynomials in /C to 
be included in massive GR. For this, we replace in (9-11) the matrix elements 11^ 
by the matrix elements K% defined in (6). As a result of this procedure, we get the 
Lagrangian density [1]: 

£ = (R + ™ 2 (4eto + ^ECC) + (/C))) . (12) 

Since all terms in (8) with n > 5 vanish identically, by construction all terms C^J r 
with n > 5 in (12) are also zero. Hence, the most general Lagrangian density (12) 
has three free parameters, to, 0:3 and 0:4. 

As it is straightforward to see, Minkowski space is a vacuum with 4> a = x a , and 
the spectrum of the theory (12) contains a graviton of mass to; the graviton also 
has additional nonlinear interactions specified by the action at hand. 



2.1 Proof of the absence of FRW cosmologies 

Let us begin by considering homogeneous and isotropic solutions to the theory (12). 
There exists a coordinate system in which the most general ansatz consistent with 
these symmetries reads as follows 

ds 2 = -dt 2 + a 2 {t)dx 2 , 0° = /(*), <p l = x\ (13) 

Here and in the following we are assuming a flat 3-d metric, but our conclusions on 
the absence of the homogeneous and isotropic solutions do not change if we allow 
for a more general maximally symmetric 3-space. Plugging these expressions for the 
metric and scalar fields into (12), and setting for simplicity 0:3 = 0:4 = 0, one obtains 
the following Lagrangian for a and /, 

C = 3M 2 ! (-ad 2 - TO 2 |/|(a 3 - a 2 ) + TO 2 (2a 3 - 3a 2 + a)) , (14) 
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where overdot denotes the time derivative 8q. We emphasize that the quantity / 
appears in the Lagrangian only linearly. The same remains true if we keep nonzero 
a 3 and a 4 - it is just the special structure of the terms £^(/C), n = 2,3,4 in (12), 
that ensures that / enters only linearly! This is a consequence of the fact that in the 
decoupling limit the equations of motion of this theory have no more than two time 
derivatives acting on the helicity-0 field in particular (and on any field in general) 
[3]. Away from the decoupling limit this is related to the constraint that was found 
in Refs. [1, 4, 5]. Here we see the constraint for the FRW metric to all orders, by 
taking variation of (14) w.r.t. /: 

m 2 d (a 3 -a 2 ) =0. (15) 

This constraint makes time evolution of the scale factor impossible. As we have noted 
above, keeping the K 3 and /C 4 terms in (12) can only modify the polynomial function 
of a on which d acts in (15). Therefore, there are no nontrivial homogeneous and 
isotropic solutions in the theory of massive GR, defined by (12). 

It is also instructive to show the absence of FRW solutions in the unitary gauge, 
for which <p a = S^x 1 *, and no / field appears in the action to begin with. In this gauge, 
the most general homogeneous and isotropic ansatz involves the lapse function N(t), 

ds 2 = -N 2 (t)dt 2 + a 2 (t)dx 2 , (16) 

and the Lagrangian (12) with a 3 = a 4 = reads 

£ = 3M 2 , {-^- - m 2 (a 3 - a 2 ) + m 2 N(2a 3 - 3a 2 + a)^j . (17) 

As can be straightforwardly verified, the condition (15) in this case arises as the 
requirement of consistency of the equations of motion for the two fields, a and N 
in (17). More specifically, one can obtain (15) by taking the difference between the 
time-derivative of the e.o.m. for N and the e.o.m. for a. Technically, this is so 
because the second term on the r.h.s. of (17) has no factors of N in it and the 
constraint arises as the direct result of the Bianchi identity of GR. 

We briefly note that the homogeneous and isotropic solutions would not be for- 
bidden if the mass term were not an explicit constant, but instead emerged as a 
VEV of some field-dependent function; i.e., if we replaced m 2 — > m 2 (a) in (12), 
where a is a scalar field that also has its own kinetic and potential terms. Then, 
variation w.r.t. / would give rise to a constraint 

d (m 2 (a)(a 3 -a)) = 0, (18) 

that relates time evolution of the scale factor to that of the a field, but it does not 
forbid homogeneous and isotropic solutions. Hence, the absence of the homogeneous 
and isotropic solutions is an intrinsic property of massive GR with an explicit mass 
term, as in (12). By this property it could potentially be distinguished observation- 
ally from the theory with a dynamical mass m 2 (a). Moreover, for the latter theory 
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one should expect the presence of an additional massless (or very light) scalar at 
cosmological distances, which is hidden by the Vainshtein mechanism at shorter 
scales. One example of this is when m 2 — > m 2 exp(a/M P1 ), for which the DL theory 
(with the kinetic term for a) reduces to a theory with two galileons coupled to the 
tensor field. 

The above- described properties of massive GR are similar to those of a peculiar 
scalar field theory (the so-called Cuscuton), defined by the following Lagrangian [21]: 



where \i is some dimensionful constant. Assuming = (pit), and the homogeneous 
and isotropic FRW metric (13), the scalar field equation reduces to a constraint, 
similar to (15) 



Therefore, the theory (19) does not possess homogeneous and isotropic cosmological 
solutions in full analogy to massive GR described above. Pursuing this analogy 
further, the homogeneous and isotropic cosmological solutions would be permitted 
if we were to promote the parameter \i into a field dependent function /i 2 — > /U 2 (o"). 
In this case the constraint would read as d (fi 2 (a)a 3 ) = 0; the latter links the time 
evolution of the scale factor to that of a, but it does not forbid homogeneous and 
isotropic cosmological solutions. 

We note that the equations of motion of the theory (19) are invariant under 
the replacement — > E(4>), where E is an arbitrary differentiable function. Hence, 
by finding a particular solution one immediately generates an infinite number of 
solutions. Whether a somewhat similar invariance exists in the equations of motion 
of massive GR is not obvious. 

Last but not least, we complement the present section with the discussion of the 
degravitating solution of massive GR (for degravitation see [22] , for the correspond- 
ing decoupling-limit solution, see [15]). Although the constraint (15) forbids the 
time-evolving FRW solutions, it still allows for a static Minkowski metric - even in 
the presence of a Cosmological Constant. To see this, we note that the action (17) 
in the presence of the vacuum energy density £ has an additional contribution of 
the form —Na 3 £. In this case one can determine the value of the scale factor from 
the iV-equation, which for a static homogeneous field reduces to, 



where, as before, p co = 3M|jm 2 . The value of N is determined from the equation 
of motion for a. The solution to this equation exists if the following inequality is 




(19) 




(20) 



2a 3 - 3a 2 + a = a 3 — , 

Pco 



(21) 



7 



satisfied, 
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(22) 



Hence, the degravitation works for an arbitrarily large positive vacuum energy den- 
sity. If we were to include also an arbitrary a 3 and a 4 , degravitation could have 
been achieved for arbitrary £. However, we note that fluctuations on the degravi- 
tated background exhibit the Vainshtein mechanism at the scale determined by the 
degravitated vacuum energy - the larger the degravitated energy, the smaller is the 
corresponding Vainshtein radius [15]. Because of this, there is an unscreened fifth 
force, and we cannot be living today on such a background. Hence, the degravitation 
could have only taken place in a far past after which the universe must have transi- 
tioned to a different background (see more in [15]). Note that the screening solution 
exists for a broad class of external sources, not just for a Cosmological Constant. 

3 Cosmology at high densities: p ^> p co 

We begin by thinking of the universe filled with pressure- less dust of density p ^> p co , 
with p co defined in (2). As was discussed in Section 1, cosmological evolution of such 
a universe can very well be approximated by the standard FRW metric of GR. This 
is so because an arbitrary Hubble patch enclosed by a sphere of radius H~ x = 
(p/3Mp l )~ 1 / 2 , is well within its Vainshtein radius, H^ 1 <C r*. Hence, the early 
universe in massive GR would evolve as it does in GR, with some small corrections. 
These corrections, for any observer in such a universe, can be estimated as some 
positive power of the ratio (m/H) <C 1. On the other hand, we would expect the 
scalar fields, a , to be in a nonperturbative (Vainshtein) regime at these scales, and 
yet, their stress-tensor should be sub-dominant to the matter /radiation stress-tensor 
that drives the FRW expansion. That this is so is a necessary condition for the self- 
consistency of the solution. Below, we will calculate the expressions for the scalar 
fields 4> a in an FRW background, and discuss its backreaction. 

The most general spherically symmetric solution, including the four scalars a , 
can always be put in the following form, 



The advantage of the latter form of the metric is that it is easier to compare to the 
standard FRW, while the appropriate 4> a fields can be treated separately. 

The Einstein equation, obtained by varying (12) with a 3 = a 4 = with respect 
to the metric, reads as follows, 



ds 2 = -dt 2 + C(t, r)dtdr + A 2 (t, r) [dr 2 + r 2 dQ 2 ] , 
<f>° = f(t,r), <P = g{t,r)-. 



(23) 



(24) 
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where TjifP is the effective stress tensor due to the mass term in (12), while T f 



(m) 



denotes the stress energy tensor of standard matter. Taking a covariant derivative 
of the above equation leads to the Bianchi constraint, m 2 V 'T^fP = 0, which is just 
the equation of motion obtained by varying the action w.r.t. <f> a , 



5S 

5<f) a 



0. 



(25) 



As discussed above, we will be neglecting at the zeroth order the m 2 T^ term, as 
well as the (possible) r dependence in the metric (23). Hence, the zeroth order 
solution for the metric will be the standard FRW solution, with A(r,t) = a(t) and 
C(t,r) = corresponding to the matter content encoded in For the scalars 

<p a , instead, we have to solve the full equations (25) in the background FRW metric 
just defined, since these are already proportional to m 2 . For this, we can rewrite 
the potential part of the massive GR action (7) using the following identity, 



KtKl - {K a a f = -12 + 6 tr + tr g^E - (tr v^S) 2 , (26) 

where the matrix ((? _1 S)^ is defined as follows: 



(p 


-g 2 


//' - gg' 





\ 


99 


-if 


-f' 2 +9' 2 










a 2 


a? 










9 2 
a 2 r 2 





V 











9 2 ) 

a 2 r' 2 / 



(27) 



Evaluating the eigenvalues of the latter matrix, and varying the resulting action 
w.r.t the two fields f(t,r) and g(t,r), one obtains the corresponding equations of 
motion, 



y ( !■ // 



— 8 r 



y i f ■ \ 2 2 ■ 

+ jig + fxa r g 



0, (28) 



dt 



y ( ' i f \ 2 2 si 



8 r 



a^fX V a 



6a 2 r + 2ag + 2a 2 rVx = , (29) 



where dot denotes time derivative and prime denotes derivative w.r.t. r, and we 
have introduced the following notation, 



g' 



y = 2ga'r-3aV, X = [ f + ^ J - Ig + ] , n = sgn(fg'-gf') . (30) 



By solving (28)-(29) for the Stiickelberg fields, and calculating the effective stress- 
energy tensor TfuP on the solution, one can evaluate the backreaction on the geom- 
etry from the presence of a 's. For the Vainshtein mechanism to be operative, two 
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conditions should be met inside the Vainshtein radius: (I) the backreaction of the 
Stiickelberg stress-tensor should be negligible, so that the background evolution is 
described by the FRW geometry to a very high precision, and (II) the metric fluctu- 
ations should be those of GR to a very high precision. This would guarantee, that 
the scalar is successfully screened. Generically one would expect these properties to 
hold in the Vainshtein region as the stress-tensor of the Stiickelberg fields in multi- 
plied by a small parameter m 2 . Hence, one should anticipate significant departures 
for the standard GR results at scales of order 1/m. Even though the metric is homo- 
geneous and isotropic, the Stiickelberg fields are not. Due to these inhomogeneous 
fields there is a physical center in each 1 / m-size domain. This center is not felt by 
matter coupled to the metric, but perturbations will be sensitive to it. This by itself 
restricts the value of m to be an order of magnitude smaller than H , or less. 

For illustrative purposes, we give in Appendix A a particular exact solution, 
which satisfies the condition (I) in an interesting way. On that solution the stress- 
tensor T^ v exactly coincides with the stress-tensor of dark energy with the energy 
density ~ p co , in spite of the fact that the Stiickelberg fields are inhomogeneous! As 
a result, this solution can exist even when no external stress-tensor is introduced. 
Hence, it is in a class of self-accelerated solutions. This solution can exist not only 
in the Vainshtein regime, but also outside of it. 

What is however not clear is whether the fluctuations on this self-accelerated 
solution are close to those of GR (it is easy to show that subhorizon fluctuations are, 
but one needs to demonstrate it for larger scales as well, which requires some careful 
calculations). Until this is known the solution can only serve a demonstrational 
purpose showing the smallness of the backreaction in the Vainshtein regime. This 
is precisely how we regard this solution in the present work. 

Putting the issues of perturbations aside, on the self-accelerated solution the 
value of the mass is related to that of the present-day Hubble parameter as m 2 = 
C Hq, where C is a free constant in the theory, which would depend on the param- 
eters o 3 and o 4 (if we were to include them); without significant tunings of these 
parameters the theoretical value of C should not be assumed to be outside of the 
interval C ~ (0.01 — 1). In this case, it should be possible to rule out such a scenario 
observationally (or at least to rule out a significant fraction of the parameter space 
for 0:3,0:4), as at the present-day Hubble scales one would expect departures from 
the FRW evolution of the order of C. 



4 Cosmology at low densities: p p co 

After the energy density p drops below its crossover value (2), massive gravity enters 
the linear regime, as discussed in Section 1. In this regime, no matter how small 
the graviton mass, the massive theory differs from the massless one by quantities 
of order 0(m ), thus exhibiting the vDVZ discontinuity [14]. Therefore, the cos- 
mology described by the massive theory is expected to differ significantly from the 
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conventional one. The purpose of this section is to study that cosmology. For this 
we first recall the status of linearized cosmology in GR. There are some subtleties 
in this, and we would like to emphasize those relevant for our discussions. 

First, the matter /radiation stress-tensor should be conserved in this approxi- 
mation, d^T^ = 0. Then, if we were to choose a diagonal stress-tensor, = 
diag(—p,p,p,p), the conservation would impose p to be time- independent. To avoid 
this restriction, we will have to choose a coordinate system in which the stress-tensor 
is not diagonal and takes the form: 

T ( P -H(p + p)x i \ , . 

^ v \-H(p + p)x l pb l3 J' [6L) 

where H, p, p are arbitrary time- dependent functions. It is straightforward to check 
then that the condition d^T^ = leads to the proper conservation equation, p + 
3H(p + p) =0, as soon as H,p and p are interpreted as the Hubble parameter, 
density and pressure respectively 3 . 

Second, strictly speaking, linearized GR itself has no homogeneous and isotropic 
cosmology: indeed, assuming that h^ v is a function of t only, the 00 component of 
the Einstein tensor vanishes, and the 00 component of the Einstein equation cannot 
be satisfied. The way out is to resort to the Fermi coordinate system in which the 
metric takes an inhomogeneous form (for recent discussions see, e.g., [24]): 

ds 2 = - (l - (H + tf 2 )x 2 ) dt 2 + (l - \H 2 ^j dx 2 = + /^ RW ) daW, (32) 

where the corrections to the above expression are suppressed by higher powers of 
iJ 2 x 2 . As long as H 2 x. 2 <C 1 , the above metric describes cosmology in any local 
patch as a small deviation from Minkowski space. 

In GR there exists a coordinate transformation that brings (32), in the ap- 
proximation considered, to a homogeneous and isotropic metric in a comoving co- 
ordinate system t c , x c , in which the stress-tensor has a conventional form, T^ v = 
di&g(p(t c ),5ija 2 (t c )p(t c )). This coordinate transformation takes the form (see [11]) 

tc = t-^H(t)x 2 , x c =^(l + ^H 2 (t)x 2 ^j , (33) 

however, the transformation itself is essentially nonlinear 4 . 

Linearized massive gravity is not much different in that respect - it does not 
admit homogeneous and isotropic solutions either. We prove this by assuming the 

3 For momentum conservation one should use the full covariant expression V^T^i that leads to 
the acceleration equation; in this case however, it is obtained at the linear order in x and requires 
correcting the stress-tensor by 0(H 2 x 2 ) quantities [23]. 

4 This is consistent with the fact that the Friedmann equation in the comoving system relates 
the square of the perturbation to density, {8a) 2 ~ GnP, where 5a denotes a departure of the scale 
factor from the Minkowski space, 5a = a — 1. 
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opposite and showing the contradiction. For this, consider the Fierz-Pauli (FP) 
Lagrangian to which any consistent Lorentz-invariant massive gravity should reduce 
at the linearized level [25] , 

C = -h^Sph pa - l -m\h^h, u - h 2 ) + W V T, V . (34) 

Here, h^ u = g^ u — r]^ u and S denotes the linearized Einstein operator 

££V = " d»d a h au - d u d a h afl + d&h - -q^Uh + v^d a dph a P) , (35) 

all indices are contracted with the flat metric, and the Planck mass has been set to 
one. Furthermore, applying & 1 to the equation of motion obtained from (34) and 
using the Bianchi identity, one gets a constraint, 

= d u h. (36) 

In the unitary gauge the metric perturbation represents a physical field, and 
requiring homogeneity makes all of its components space- independent, 

di V = 0. (37) 

With this assumption, the constraint (36) reduces to the following equations 

hoi = const = Ci, ha = const = c. (38) 

The second of these equations implies that either a solution is trivial or else it 
cannot be isotropic. Hence, linearized massive GR has no nontrivial homogeneous 
and isotropic solutions. Moreover, it is straightforward to check that in general, 
inhomogeneity cannot be removed entirely into the longitudinal degrees of freedom. 

Then, what is a snapshot of the universe at p p co ? We will show below that 
it can be pictured as a collection of multiple domains, each of size 1/m, such that 
well within a given domain, at scales < 1/m, cosmology deviates significantly from 
the conventional GR one. How these domains are glued together is a complicated 
question that is not addressed here. Nevertheless, at a scale much greater than 1/m, 
when averaged over many domains enclosed by this scale, the universe should look 
homogeneous and isotropic again. 

What exactly is then a solution to the linearized massive GR in each domain? 
To get this solution, it is useful to introduce appropriately normalized Stiickelberg 
fields: 

V = V + %v v + d u v, = v + d " Av + dvA > + > ( 39 ) 

m m z 

and consider scales that are much smaller than 1/m. In this approximation, the FP 
theory (excluding the totally decoupled vector mode) reduces to: 

C = -h^S^h pa - h^(d,d u n - t^Dtt) + h^T, u + 0(m 2 ) . (40) 



12 



The latter can be diagonalized by the conformal shift h^ v = + t/^tt, giving rise 
to the Lagrangian 

£ = -\w v £ p °h p<J + ^TrDvr + hTT^ + vrT + 0(m 2 ) . (41) 

The equations of motion that follow from this Lagrangian are: the GR equations 
for h^ u , and a simple equation for n, Dn = — T/3. The solution for h^ u , as argued 
above, is h F J} w given in (32), while for ir we obtain n-K so i = 2(H + 2H 2 ). As a result, 
the physical metric is h P ^ ys = h F J} w + r]^ u Tc so i, and the interval takes the form: 

ds 2 = -(l-±(2H + # 2 )x 2 ) dt 2 + (\ + i(2if + H 2 )^j dx 2 . (42) 

Note that the linearized Ricci curvature on the physical metric h P ^ ys is zero (in the 
leading approximation), as it should be the case for all cosmologies due to (36). This 
metric deviates from its GR counterpart by the value of 7r so i, which is of the same 
order as h FPW itself. This is a cosmological manifestation of the vDVZ discontinuity. 

At scales ~ 1/m the mass terms neglected in (41) should be reinstated and used. 
However, at yet larger scales, which enclose a large number of domains, the universe 
should look homogeneous, if we average over all enclosed domains. We notice that 
at scales ^ 1/m, all the derivative terms in the Lagrangian (34) should be neglected 
and only the mass terms should be kept. Then, the equation of motion takes the 
form: 

m 2 ( V - n^h) = 2(7^) , (43) 

where (T^) is the stress-tensor (31) averaged over many 1/m-size domains. At such 
large scales gravity is screened. Depending on the initial conditions there may be 
a number of different scenaria of how one could match the large and small scale 
behavior to each other. If the 1/m-size domains are densely packed and adjacent 
to each other, then the matter/radiation will get accumulated near the boundaries 
of the domains, as it is expelled from the bulk. However, if the domains are well 
separated, or there is only one domain, then density near the boundaries should be 
expected to be suppressed since gravity of the bulk material is screened and matter 
particles will free stream out of the domain. All these scenaria, and the initial 
conditions that could give rise to them, need separate detailed studies. 

Acknowledgments 

We would like to thank Matt Kleban and Roman Scoccimarro for useful comments. 
The work of GD'A is supported by the James Arthur Fellowship. CdR is supported 
by the Swiss National Science Foundation. The work of SD is supported in part 
by the NSF grant PHY-1068438. SD thanks the hospitality of the Kavli Institute 
for Theoretical Physics at UC Santa Barbara (under the NSF grant PHY05-51164), 



13 



whilst this work was being completed. GG is supported by the NSF grant PHY- 
0758032. The work of DP is supported by Mark Leslie Graduate Assistantship. AJT 
would like to thank the Universite de Geneve for hospitality whilst this work was 
being completed. 



A An exact background solution 

We will now show that it is possible to find a solution of the system (28)- (29) 



a 



y i , 9 \ , 22/ 
T V ~a ' ^ r 9 



+ fig + fxa r g 



0, (I) 



ay/X \ a J 

- 6a 2 r + 2ag + 2a 2 rVx = . (II) 

We notice that (I) is identically satisfied for y = 0, or 

g(t,r) = ^a(t)r. (Ill) 
For this ansatz eq. (II) reduces to the following equation for f(t,r): 

\[{^+f) -(l Har + fI ^) =^' + l~ » Hrf ' ' (IV) 

with H — a/a denoting the usual Hubble parameter. 

The structure of the equation (IV) guarantees that all square roots appearing 
in it are always well-defined. Indeed, the above equation for / can be reduced to a 
simpler one by rewriting it as follows, 

^{^f + \- ^ Hr f) ~ \HW - S + 2Hrff> - = 

pf + 1 - nHrf. (V) 
Squaring the latter equation puts it in the following form: 

C (1 + a 2 H 2 r 2 ) - 2Hrff + -a 2 H 2 r 2 = . (VI) 
a 2 4 

Obviously, the expression under the square root in (V) is positive semi-definite for 
any solution of the squared equation. Moreover, for any solution of this equation, 



t 



y_ 
x 



g + H—)+Liarf 
a 



— 8 r 
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the quantity appearing on the r.h.s of (V) is positive, which can be seen from writing 
it as 

Pt +1- *Hrf jf) + 3 -, 

and recalling the definition of /i, eq. (30). Therefore, any solution of eq. (VI) will 
solve eq. (IV). 

One solution of eq. (VI) is 



f(t,r) = 



9 



(VII) 



16T J a(t)H(t) 

where T is an integration constant with dimensions of time, and the choice of the 
lower limit of integration corresponds to a constant shift of 0°. 

We are now in a position to compute the stress-tensor of the 4> a fields and compare 
it with that of matter. It follows from (7), that the effective stress-tensor is given 
by the following expression: 

1 6 W=g 



T (K) 



-9 8g 



(l2 + (Tr - 6)Tr - <? Q/3 E Q/3 ) 9ilv 
+ 2£^ + (3 - Trv/^FS) (g^y/rWv + du^^F^Y^ 



The nonzero components of Tf£P are therefore: 



■00 



12- I \ \ | — - - 

ar 2 



ar \ar / 



3 -'-H/r f 2 '-(!<>' !,.f 

X \ ox J \ a 



(VIII) 



and 



Ti5 } = -a 2 



12 + 4VX( 2 A (A_ 6 

ar 2 1 ar \ar 



('-'£)((*)"-(£)"■►*(>'-'') 



(K) 



12 + 4Vx(°-- 3 -)+2°- fA_ 6 
\ar 2 J ar \ar 

a V / \ar / ar V ar/ 



sin 2 9 2 
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Remarkably enough, as we show below, for the solution at hand the inhomogeneities 
of the Stuckelberg fields completely fall out from the expression for the effective 
stress-tensor T$\ 

Using the exact solution for g(r,t), as well as the equation of motion for f(r,t) 
(VI), exactly reduces to the diagonal, cosmological-constant-type form - with 
the corresponding Hubble scale set by the value of the graviton mass, 



T {K) 



V 



/ 


-p 








0\ 







p 
















P 





V 











PJ 



P = -P= \™ 2 Ml x = l -p co . (IX) 



Therefore, the backreaction from the Stuckelberg fields is indeed negligible for the 
universe filled with matter or radiation with density significantly exceeding ~ p co . 

The background solution we just found, consisting of the FRW metric plus eqs. 
(Ill) and (VII), is exact. Cosmology for this solution - at least at the background 
level - is therefore completely insensitive to the presence of inhomogeneity in the 
Stuckelberg scalars. Even in the absence of any external sources, the geometry 
describes the homogeneous and isotropic self-acceleration of the universe with the 
Hubble constant equal to m/2, all of the inhomogeneities of the solution being 
removed into the Stuckelberg scalars. On the other hand, the space-dependence of 
the background value of <p a can be probed by perturbations on the FRW metric. 

More generally, we expect the theory to admit solutions with truly inhomoge- 
neous geometry - with the metric being impossible to put in a homogeneous form 
by any coordinate transformations. For such solutions, we expect the backreaction 
of the Stuckelberg fields to be negligible in the high-density regime, while becoming 
important for densities ^ p co . We emphasize again that the question of whether 
the evolution can continue on the self-accelerated solution, or alternatively should 
switch to the low density vDVZ regime, remains open, as the matching between 
these two regimes is hard to study analytically. 



B Anisotropic solutions 

In this Appendix, we show that there are cosmological solutions which are homoge- 
neous but anisotropic. 

Parity symmetric solutions 

Let us start by considering the following diagonal Ansatz for the unitary gauge 
metric: 

ds 2 = -N 2 (r)dr 2 + a 2 (r)dx 2 + a 2 2 (r)dy 2 + a 2 3 (r)dz 2 , (X) 

which is the most general homogeneous metric invariant under the discrete parity 
symmetry x — > —x. Before plugging this into the action, it is convenient to redefine 
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the time variable as 

dt = A(r)dr . (XI) 

This gauge transformation will excite one of the Stiickelberg fields, which will now 
read: 

0° = /(f) , f = x\ (XII) 

In this gauge we can immediately derive a constraint on the scale factors. 
The matrix = g^ a d a (j) a d u (j) a is diagonal, with the following eigenvalues: 

111 
o"o = / , fi = — , 02 = — , 03 = — • (XIII) 

(X-^ 6^3 

The mass term in the action therefore reads 

— 12a!a 2 a3 + 6(aia 2 + aia 3 + a 2 a 3 ) — 2(ai + a 2 + 03) 



5 = Mk m 2 I d 4 x 



+ 2/(3aia 2 a3 — aia 2 — 0103 — a 2 as) 



• (XIV) 



It is immediate to derive the following algebraic constraint on the scale factors, 
which is proportional to m 2 and it is valid independently of the matter Lagrangian 
coupled to the metric: 

3aia 2 a 3 - (aia 2 + aia 3 + a 2 a 3 ) = k , (XV) 

where k is an integration constant. As a check of the calculations, we can easily see 
that there is no isotropic solution other than Minkowski space-time; for a x = a 2 = 
a 3 = a, we have a 3 — a 2 = k, which is exactly eq. (15). 

It is simple to show that we cannot have a solution for which all the scale factors 
grow at small times. In fact, requiring that a* — > as t — > 0, we see that we should 
set k = 0, because the whole LHS should vanish for small times. However, since the 
scale factors are positive, the equation cannot be satisfied. Thus, we cannot have a 
solution with all the scale factors growing with time. The only non-trivial solutions 
will have one direction contracting and the others expanding, or vice-versa. 

Axisymmetric solutions 

We now discuss the most general homogeneous and axisymmetric solution. The 
most general anisotropic ansatz can be written as: 

<i>° = f(t) + b^, (XVI) 
0* = A\x j + c\t) , (XVII) 
ds 2 = -dt 2 + a 2 (£)da; 2 + a 2 2 (t)dy 2 + a 2 3 (t)dz 2 , (XVIII) 

where bj is a constant vector and A 1 - is a constant matrix. 
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If we impose axial symmetry around the third direction, our ansatz should satisfy 



k = (0, 0, q) , A\ = diag(l, 1, B) , c\t) = (0, 0, c(t)) (XIX) 
ai(t) = a 2 (t) = a(t) , a 3 (t) = b(t) . (XX) 

The matrix = g aX d\(f) a d u (f) b ri ab is block diagonal and can be easily diagonalized. 
Its eigenvalues are given by: 



1 2 9 + vV ~ h 2 9- Vd 2 ~h 

"I "2 

where 



2 2 a ' V » ^.2 » V » /VYI^ 
a 2 = "2 ' °3 = 1 » ^4 = 7 , (XXI; 



1 

9 



^-q 2 + b 2 {f 2 -?)\ , h = b 2 (qc-Bf) 2 . (XXII) 

We can derive two constraints by varying the action with respect to / and c. After 
some algebraic manipulations, we arrive at 



' B + q-fjbjf + c) = a 2 (B + q)-^ 

B - q - fib(f - c) b(3a 2 -2a) ' 1 ' 
,2 [a 2 (B-q)-fik + ][a 2 (B + q)-fik^ 

b = (3a* - 2af ' (XXIV) 

where \i = sgn(gc — Bf). 

If we require that both a and b vanish when t — > 0, we must choose k+ = k- = 0. 
Now, from the second constraint, we find B 2 > q 2 . In order to have a positive 6, 
the only possible choice would beB + q<0^>B<0. We can prove that this 
is an inconsistent solution by looking at the square root. First, we notice that the 
symmetry <p a — > —<p a allows us to study only the case in which B + q — (j,b(f + c) > 0. 
So, we would like to have 

fib(f + 6)<B + q<0, nb{f -6)<B-q<0. (XXV) 

Now, these conditions have no solution. In fact, they imply 

nf < , (XXVI) 

which is absurd. Indeed, we also see from (XXV) that |/| > |c|, which along with 
\B\ > \q\ and B < implies that \i = sgn/, so necessarily \ij > 0. In conclusion, 
also in this case we cannot describe a Universe which expands starting from a small 
initial volume at early times. 
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